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Abstract 

In the semi-classical limit the relativistic quantum mechanics of a stationary beam 
of counter-streaming (negatively charged) electrons and one species of positively 
charged ions is described by a nonlinear system of finite-temperature Thomas-Fermi 
equations. In the high temperature / low density limit these Thomas-Fermi equa- 
tions reduce to the (semi-)conformal system of Bennett equations discussed earlier 
by Lebowitz and the author. With the help of a sharp isoperimetric inequality it is 
shown that any hypothetical particle density function which is not radially symmet- 
ric about and decreasing away from the beam's axis would violate the virial theorem. 
Hence, all beams have the symmetry of the circular cylinder. 
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1 INTRODUCTION 



Modern books on charged-particle beams, e.g. [^T|], usually contain a chapter about the 
Bennett model [|J, but back in the early 50's when regular research on charged-particle 
beams came into sharper focus, W.H. Bennett's pioneering pre- WWII paper j|] on the 
statistical mechanics of a relativistic, stationary particle beam had been forgotten, appar- 
ently, and so in 1953 Bennett sent out a reminder note [H. For some reason or other, 
Bennett's note did not appear until 1955 ||, the very year when Joel L. Lebowitz was 
launching his stellar career pOf with center of gravity in stationary non-equilibrium statis- 
tical mechanics [31], |32], |33|. At that time, a single issue of The Physical Review was still 
of a decent size and could be consumed from first to last page by an individual with huge 
scientific appetite such as Joel, and Bennett's note || did not pass unnoticed before Joel's 
hungry eyes. All this happened a few years before I was born, but when I came to spend 
some postdoctoral time with Joel nearly 40 years later, several interesting questions raised 



by Bennett's work were still unanswered, and so we began to answer some of these |29 
One of the problems we had to leave open was that of the symmetry of a beam. Following 
Bennett we only inquired into circular-cylindrically symmetric solutions. While it is a nat- 
ural conjecture that in the absence of external fields an unbounded straight particle beam 
with finite electrical current through its cross-section necessarily posesses the symmetry 
of the circular cylinder, how to prove it is not quite so obvious. It is with great pleasure 
when in this paper I present a rigorous proof to Joel. 

Fitting for the occasion, the proof of the cylindrical symmetry of the beam involves 
statistical mechanics in an essential way. Namely, it is shown that any hypothetical sta- 
tionary beam with finite electrical current whose particle density functions are not radially 
symmetric about and decreasing away from the beam's axis would violate the virial theo- 
rem for this many-particle system. This symmetry proof covers Bennett's strictly classical 
model as well as its semi-classical upgrade, i.e. a system of relativistic, finite-temperature 
Thomas-Fermi equations which in the high-temperature / low-density limit reduce to the 
(semi-)conformal Bennett equations. The proof is, however, restricted to a system of 
merely two equations because the coefficient matrix for the beam equations has rank 2. 
Our symmetry theorem therefore does not apply to beams that consist of the negatively 
charged electrons and more than one, differently positively charged ion species; but then 
again, our method of proof not only yields the cylindrical symmetry of the beam, it also 
yields monotonic radial decrease of the particle density functions. Hence, it is conceiv- 
able that monotonicity of the density functions may be violated in an electron / multi-ion 
species beam while cylindrical symmetry might still hold — yet to prove that would seem 
to require an entirely new argument. 
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Incidentally, our result also sheds some new light on the theory of white dwarfs ||. 
These Earth-sized, expired stellar objects shine in bright white light because they are still 
incredibly hot compared to our Sun, yet they are relatively cold compared to their Fermi 
temperature and therefore essentially in their quantum ground state. This justifies using 
zero-temperature Thomas-Fermi theory for the description of their overall structure H] - 
a fortunate happening, for finite-temperature Thomas-Fermi theory could not be used in 
three dimensions since it does not have solutions with finite mass. Interestingly, the finite- 
temperature Thomas-Fermi equations of the two-dimensional caricature of such a white 
dwarf star should have solutions with finite mass, because the gravitational potential in 
two dimensions is sufficiently strongly confining for this purpose. In any event, relatively 
little is known rigorously^] for such a gravitating plasma of negative electrons and positive 
nuclei (all species treated as fermions) in either two or three dimensions; see the discussion 



of this model by W.E. Thirring in the preface to the E.H. Lieb jubilee volume [pflj| , where 
Thirring gives an amusing account of the pitfalls associated with the fact that the Thomas- 
Fermi equations are the Euler-Lagrange equations for the saddle points of a variational 
functional. When dealing with saddle points, existence and symmetry of solutions via 
minimization by radial decreasing rearrangement [[I], || [7|, ^] is not an option, and neither 
is symmetry via uniqueness by convexity [^7j] of the functional. 



Now recall that by the Biot-Savart law the magnetic interactions of straight, paral- 
lel electrical current filaments are attractive, with a distance law that is identical to the 
Newtonian gravity law in two dimensions. From this it follows that the finite-temperature 
Thomas-Fermi beam equations are identical to the finite-temperature Thomas-Fermi equa- 
tions of the two-dimensional caricature of a white dwarf model, with the magnetic flux func- 
tion re-interpreted as the gravitational Newton potential in two dimensions, and the mean 
electric current of each species (positive after at most a joint space rotation) re-interpreted 
as the mass of that species. Our symmetry result can be rephrased thus: two-dimensional 
finite-temperature white dwarfs are radially symmetric. 

Our proof of symmetry, which is based on the Rellich p9| -Pokhozaev p8| identity 
(which expresses the virial theorem) and the classical isoperimetric inequality [[|, [|, does 
involve radial rearrangements in a strategy that goes back at least as far as H, where it 



^^More is known rigorously [|[ |35| for the locally neutral approximation of the model, where the positive 
and negative charges are distributed identically and Coulomb's law is discarded. In particular, radial 
symmetry of solutions for this locally neutral model has been proven by energy minimization through 
radial rearrangement ]35| ] . We remark that due to the enormous ratio of the electrical and gravitational 
coupling constants the locally neutral approximation is expected to be an excellent approximation for a 
white dwarf; however, this is not generally the case for a particle beam, where the ratio of electric and 
magnetic coupling constants may be arbitrarily close to 1. 
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is applied to Liouville's equation?] in a disk cR 2 0. In |IJ this strate gy was generalized 
to systems of PDEs of Liouville type in all IR 2 which are unrestricted in size but which 
have a symmetric, fully stochastic coefficient matrix of full rank. The Bennett equations 
also constitute a Liouville system, but are not covered by the theorem of because their 
coefficient matrix is generally not symmetric, has some negative elements, and is always 
rank 2. The present paper develops the necessary generalizations of ]10 ] to overcome the 
first two peculiarities of the Bennett equations, but the rank 2 restricts the proof to a 
system of two equations. By adapting the treatment of single PDEs with more general 

]) and Jy], 0] to the system case we are 



nonlinearities developed in |36 



(cf. also 

able to extend our proof of symmetry for the Bennett equations to the relativistic, finite- 
temperature Thomas-Fermi beam equations. 

Our proof simplifies considerably when the systems of Thomas-Fermi and Bennett 
equations are restricted to a disk with O-Dirichlet boundary conditions for the electric and 
magnetic potentials. In this compact case, an alternate proof of the radial symmetry and 
decrease of the solutions to systems of PDE which includes the finite-temperature Thomas- 
Fermi and Bennett equations, was given by Troy [fffil , who exploited Alexandroff's method 
of moving planes. For more on the moving-planes method, see jll], [19], |34|, [T3|, || |i~3|j . 
Troy's proof has been extended to Liouville systems in unbounded domains, the Bennett 
equations not included though, 



in 
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Presumably, the moving planes method can be 
made to work also for the system of Thomas-Fermi equations studied here; however, this 
is not done in this paper. 

While the present paper addresses only the question whether invariance of the PDEs 
under rotations implies radial symmetry of their solutions, these PDEs feature other sym- 
metries which deserve mentioning. The system of Thomas-Fermi equations is invariant 
under the isometries of Euclidean space, simple gauges, and Lorentz boosts along the beam. 
The Bennett equations are in addition to that invariant under isotropic scaling in M. 2 , and 
for a special family of parameter values also under Kelvin transformations, in which case 
they are invariant under the Euclidean conformal group of M. 2 . In this fully conformal case 



the conformal orbit of the finite current solutions is connected and itself invariant 29 



Invariance under the Euclidean conformal group holds also for the Liouville systems stud- 
ied in [plj, but their conformal orbit of finite mass solutions is generally not connected, 
and each component not invariant under inversions jT5]. Toda systems in M 2 , which are 
Liouville systems with symmetric coefficient matrix given by the SU(N) Cartan matrix, 
are studied in [22, 23]. The distribution of negative and positive signs in the SU(2) Cartan 

2 The elliptic Liouville equation, known from two-dimensional differential geometry, is meant and not 
the evolution equation on phase space known from statistical mechanics. 
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matrix is opposite to that in our Bennett equations, and sure enough, our radial symmetry 
proof fails in this case. Interestingly, in this case one can show that radial symmetry is in 
fact broken by some solutions, see the bifurcation argument with n = 2 in (1.7) of [ID], and 
see [ 23 1 for the construction of the complete solution family with finite masses. Another 



interesting topic not discussed further here is whether the translation invariance along the 



beam can be broken, as is suggested by various dynamical beam instabilities [44 



The remainder of this paper is structured as follows. In the next section we formulate 
the basic equations of the semi-classical beam model and its classical limit. Existence of 
solutions is briefly touched upon. In section 3 we state our two main theorems, and in 
section 4 we present their proofs. 



2 RELATIVISTIC BEAM EQUATIONS 

We let a G S 2 denote the fixed axis of the beam, x G M 2 a point in the cross-section of the 
beam containing the coordinate origin, and p G M 3 the kinematical particle momentum. 
The self-consistent electric field of the beam is given by E(x) = — V0(x), where <fi is the 
electric potential, and the magnetic field by B(x) = Vip(x) A a, where ip is the magnetic 
flux function. The beam consists of spin 1/2 electrons (negatively charged, thus indexed 
by s = — ) and one species of positively charged spin 1/2 fermions (indexed by s = +), 
characterized by the following parameters: the particle charges q s and rest masses m s ; the 
rest frame temperatures T s ; the external chemical potentials and lab frame drift speeds 
cu s , where c is the speed of light and v s G (—1,1). We demand u + ^ z/_, as appropriate for 
counter-streaming particle species. The temperatures and drift speeds combine into the 
thermal lab frame parameters = k&T s y/l — r/^. 

2.1 The semi-classical model (Thomas— Fermi theory) 

The finite-temperature Thomas-Fermi model of a straight, relativistic beam is set up as 
follows. In the lab frame the density of s-particles at x is given by p s (x) = Qj F ((p, ip) i x ), 
where 

^^) = TtJ 1 7-^- V C 1 ) 

JR 3 l _|_ e -l3s^s-cy/m 2 s c 2 +\p\ 2 +i> s cp-a-q s l<j>-i> 3 ip]) 

is the finite-temperature Thomas-Fermi density function for the relativistic s-species, 
which is subjected to the integrability condition 

/ gJ F (4>,ij)(x)dx = N s , (2) 
Jr 2 
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where N s is the number of s-particles per unit length of beam. The phase-space density 
function under the integral in ([I]) is the drifting Fermi-Dirac-Jiittner function [26] with 



local chemical self-potential —q s ((j)(x) — v s il){x)). The electric charge and current densities 
in the Poisson equations for the electric potential and the magnetic flux function ip are 
computed with the density functions (|I]), which leads to the system of nonlinear PDEs 

-A0 =^Y JS q s G7^) (3) 
-A^ = 47rE s ^J F (0,^). (4) 

Here and in the following, or ^ 4 always stands for summation over the particle species, 

i.e. s = =F and t = =p. 

The Thomas-Fermi equations (0), (f|), are invariant under the isometries of three- 
dimensional Euclidean space, Lorentz boosts along the beam's axis a, and the gauge trans- 
formation 

(j){x) -> (j){x) + O ; ip(x) -> ip{x) +^ ; fJ's ->• A* s + g s (0o -v a ij) Q ), (5) 

where <po and ipo are arbitrary constants. 

Since we are interested in the beam's natural symmetries, we will not allow "sources 
at infinity" which would deform the beam; hence, we supplement (||]) and (f|) with the 
asymptotic conditions that, uniformly as \x\ — > oo, 

lim =2= lim (6) 

|aj|— >oo Q In t— r x|-*oo / In T—r 

with I and Q ^ 0, where 7 = N s q s v s c is the total electrical current through the 
beam's cross-section and Q = J2 S N s q s the total charge per unit length of beam in the 
lab system; if Q = 0, the left equation in (||) is to be replaced by the condition that 
4>{x) — > const uniformly as \x\ — > oo. The situation I = is not considered here, for then 
of course there is no stationary beam. 

Remark: There are good reasons to conjecture that the asymptotic conditions (^|) are in fact 
implied by (|l|)-(§). Analogous results have been proven for Liouville's equation [14] and for 



some Liouville systems [I0[ [15]. No attempt will be made here to generalize these results 
to ([[])- (f|). However, we note that such a generalization would have the interesting physical 
implication (within the limits of applicability of the model) that one cannot maintain a 
stationary straight beam of finite current, whatever the geometry of its cross section, when 
there are magnetic or electric multipole sources "at infinity." 

To the best of the author's knowledge, the existence of beam solutions in the Thomas- 
Fermi model (|l])-(§) with asymptotics (§) has not yet been studied rigorously. However, 
this semi-classical model is surely more regular than the classical one, addressed next. 
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2.2 The classical limit (Bennett theory) 

In the high-temperature / low-density limit, i.e. formally < ft < 1 and f3 s ii s <C —1, the 
Fermi-Dirac-Jiittner functions |26j reduce to the Maxwell-Boltzmann-Jiittner functions 
[^4|] (see also [17], p. 46, eq. (24)), so that the Thomas-Fermi densities ([!]) simplify to 



Boltzmann densities, 
and (0) becomes 



2 

7^ 



-/3 s (c^/ m? s c 2 + \p\- 



-u s cp-a 



dpe 



I3 3 (p 3 -q s [4>-v a ip]) 



gf(M) (x)dx 



n s . 



The system of equations (|) and 

-Acj) = AnZ s N s q, 



then reduces to the Bennett equations 

e -P s qs(<t>-v s il>) 



-AtP = AirJ^sNsqsiSs 



f R 2 e-^i^-^dx 

„-j3 s q s {<t>-v s i[>) 



J R 2 e-^i^-^Hx'' 



(7) 
(8) 

(9) 
(10) 



see ||] eq.'s(8),(9), and eq.(7),[| where we have eliminated the external chemical poten- 
tials n s via (H). 

The Bennett system is invariant under the isometries of three-dimensional Euclidean 
space and under Lorentz boosts along the beam's axis, a. Restricted to the beam's cross- 
section, it is also invariant under isotropic scaling, and in the special case when the pa- 
rameters satisfy 

PMv s c- 1 I-Q) = 2, s = T, (11) 

also invariant under translated inversions. Thus, (|TTD implies invariance of the Bennett sys- 
tem under the conformal group of two-dimensional Euclidean space, acting in the beam's 
cross-section. In addition, the Bennett equations are invariant under a gauge transforma- 
tion (p(x) — > cf)(x) + cf) , ip(x) — > ip(x) + ipQ. Recall that we already eliminated the external 
chemical potentials via the constraint equations (0) in the Bennett limit. 



3 In his papers g, Bennett employed a classical, semi-relativistic setup, assuming drifting Maxwell- 
Boltzmann distributions with relativistic drift speeds, yet with non-relativistic velocity dispersion in the 
cross-section of the beam; the relativistic model with drifting Juttner functions was used in || . It should 
be noticed, though, that after integration over momentum space the very system of equations (^), ( fTo| ) 
results in either case, and it does so also in the strictly non-relativistic limit |2^| — except for minor 
re-interpretations of the parameters in each case. 
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In the conformally invariant case (pTTj) , Bennett's AnsatzQ 



x 



r x cij){x) = v{x) = q- 1 
maps (0) and ( |TU| ) separately into Liouville's equation [^7 

- Av = A-n 



L 2 e 2v dx 



(12) 



(13) 



As remarked above, it has been proven in |14[] that any regular solution of (|13"D, with the 
understanding that J exp(2t> )dx < oo, satisfies 



lim 

|a;|— »oo 



V{X) 



(14) 



uniformly as |x| — > oo, which implies that the asymptotic conditions (|]) are automatically 
satisfied if qb and ip are given by (0). It has also been proven in [H|, and subsequently in 
[p70| , pT6|| by using alternate techniques, that ( |13"D has only one regular family of solutions, 
given by 

v(x\x ;k) =v + \n -, (15) 

1 + k 2 \x — Xq\ a 

where k~ x > is an arbitrary scale length, xo the arbitrary center of rotational symmetry 
of the solution, and vq an arbitrary gauge constant. The corresponding current density 
j(x) and charge density q(x) are given by 



I- 1 



3{x) 



k 2 



7T (1 + k 2 \x - X \ 2 Y 



Q-'qix) 



(16) 



The density profile ( |T5D is the celebrated Bennett beam profile. 

Bennett speculated about the existence of other solutions to (||) and ([T0|) with asymp- 
totics (Q); see Q| p. 893, and p. 1587. (In the punctured plane additional solutions are 
readily found, see e.g. 0; however, they all lack regularity, due to a point source, at the 
origin.) In |H| we proved that in the conformal case (|TT1), Bennett's system of equations 
(|9p and fllCf), supplemented by the asymptotic conditions (||), are in fact equivalent to 
( |13"D (with asymptotic condition (0) automatically satisfied, see above) so that ( |15D then 



4 Bcnnctt actually made the Ansatz that p+{x) / 'p-(x) = const, which up to gauge freedom for the 
potentials is equivalent to (fl2l). 
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exhausts all possibilities. Moreover, for the semi-conformal case where ( |TTD does not hold 
we proved the existence of a continuous parameter family of smooth radial solutions to @ 
and (|Tt]) with asymptotics (^) which are not invariant under inversions. 

All the solutions of our beam equations are automatically also stationary solutions of 
the equations of Vlasov's relativistic kinetic theory j|3| . In [^J we showed that the Bennett 
equations can also be realized as transversal part of stationary dissipative kinetic equations 
in which the dissipation, modeled by a thermostat, compensates the action of an applied 
longitudinal electromotive force that drives the current. 

In p9| we also gave a rigorous proof that all radial solutions of (|) and fllPl) satisfying 
(H) also satisfy the Bennett identity 

c -2j2 _ Q 2 = 2J2 s N s k B T sV ^T^. (17) 

The identity (|17|) was originally obtained by Bennett || in a formal (and not entirely com- 
pelling) manner by studying the radial time-dependent virial. In this paper we will show 
that the Bennett identity ([H])> respectively it's counterpart for the Thomas-Fermi model, 
holds a priori without assuming symmetry, and this fact will be one major ingredient in 
our proof of the cylindrical symmetry of the beams. 



3 MAIN RESULTS 

To state our virial theorem, we introduce the thermodynamic potentials (per unit length 
of beam), given by 



^ TF =Eff 'I Lb* 




l _|_ e Ps(vs-cy/m%c 2 +\p\ 2 +v s cp-a-q s [tj>-v s ip]) 



dp dx (18) 



for the semi-classical model, respectively by 

jB = y^-l2 f I ^-cy/m^+W+Uscp-a-^-Us^]) dp dx (1Q) 
S h A J R 2 J R 3 

for the classical model. 

Theorem 3.1: (Virial identity.) Let <p e C 2 ' a (M. 2 ) and G C 2 ' Q (R 2 ) solve and 
(Op under the constraints (jl) ; respectively solve (j^) and (jJU) under the constraints 
s = =F, in either case subjected to the asymptotic conditions (|6[). Then 

c~ 2 I 2 -Q 2 = 2 J, (20) 

where J stands for either J TF or J B . 



9 



We also show that deviations from cylindrical symmetry violate (|20|), which is expressed 
in the next theorem. 

Theorem 3.2: (Cylindrical symmetry.) Let G C 2 > a (R 2 ) and ip G C 2 ' a (R 2 ) solve fj^) 
and (fjj under the constraints (j^j, respectively solve (^) and under the constraints 
s = =F ; subjected to the asymptotic conditions ^). Then there exists a point x Q G M? such 
that both and ip are radially symmetric about Xq, and the density functions ^(0,-0) (x) 
are decreasing away from xq, where Q s here stands for either the Thomas-Fermi or the 
Boltzmann density function. 

4 PROOFS 

We rewrite the Thomas-Fermi, respectively Bennett system in two equivalent versions, 
which may be called the "density potential representation" and the "chemical self-potential 
representation." We will switch between these representations at our convenience to obtain 
the asymptotic estimates, as \x\ — + oo, and the isoperimetric estimates needed for our 
proofs of Theorems 3.1 and 3.2. 

4.1 The alternate PDE representations 

The chemical self-potentials U s (x), x G M 2 , are given by 

U t = -q t {<f>-v a ^). (21) 

We also introduce density potentials u s (x), i6l 2 , defined by the invertible linear system 

= E s QsU s , (22) 
i 1 = Y.s^sQsUs. (23) 

Clearly, 

Us = J2t1s,tu u (24) 

where 

j s ,t = -q s qt(i ~ v s vt) (25) 
denotes the entries of the matrix of coupling constants. Notice that 

de t( T ) = -(g + g_)>+ - z/_) 2 , (26) 

so that for u + ^ z/_, we have 

rank (7) = 2, (27) 
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hence 

«. = Et%JU t (28) 

where j~ t denotes the entries of the inverse matrix 7 _1 to 7. 

Now let Q s stand for either Qj F or Qf. We note that Q s {4>ii>) depends on and ip 
only through the combination —q s {4> — v s ip) = U s ; thus we can write Q s ((j),tp) = G S (U S ) = 
G s (^2 t ls,tUt), where of course G s stands for either Gj F or Gf. In either case, the map 
w ^ G s (w) is monotonic increasing. 

It then follows at once that the chemical self-potentials U s solve the system of nonlinear 
PDEs 

-AU s = 4TrZ t ls,tG t (U t ), (29) 
supplemented by the integrability conditions 

/ G s (U s )dx = N s (30) 
and by the asymptotic conditions that, uniformly as |x| — > 00, 

lim ^ = 2^7.,^- (31) 

\x — >oo In i — r 

Alternately, in terms of the u s we get the following representation for our Thomas- 
Fermi / Bennett models, 

- Au s = 4irG s (Zt7s,tU t ), (32) 
supplemented by the integrability conditions 

/ G 8 (Ztl8,tU t ) dx = N s (33) 

and by the asymptotic conditions that, uniformly as \x\ — > 00, 

lim ^# = 2iV s , (34) 

1*1— >oo In -r-i 

\x\ 

for s = =F. This constitutes the density potential representation of our Thomas-Fermi / 
Bennett models. 

Remark: For the sake of completeness, we also state the PDEs of the Bennett model 
explicitly as a Liouville system. We readily eliminate the fx s in terms of the N s , using fl3~3"|). 
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Setting now u s = N s v s and (3 s j Sj tN t = K s>t , and furthermore ^2 t K s ,tVt = 2V, (equivalently, 
j3 s U s = 2V a ), with s and t taking the "values" ±, we rewrite fl3"2|) into the form 

exp (Y,t K s,t v t) 



Av, = 4:71 



and 



into 



AV 8 = 2tt£> 



exp (2V t ) 



s,t' 



(35) 



(36) 



'/ H a exp (2Vt)dx' 

Equations (^) and ( |36"D are explicit alternate representations of the Liouville system as- 
sociated to the Bennett model. The coefficient matrix k is manifestly non-symmetric in 
general, having negative diagonal and positive off-diagonal elements. Note that in the con- 
formal case (0), viz. J2t K s,t = 2 for s = ±, the Ansatz v + = u_ = v in (|35|) , respectively 
V + = V- = v in (|36"D, reduces both (|35| ) and (|36|) to Liouville's equation (|T3|) . 



4.2 Isoperimetric estimates 

Let G s continue to stand for either Gj F or Gf. We introduce g> s , the primitive of G s , i.e., 
g"s(w) = G s {w) for tost, such that the integrals 



/ g s {U s )dx = M s 
Jr 2 



(37) 



exist (notice that M s is defined by (|5?D). In each case this primitive g s is unique and given 
by 



gJ F (u s ) = (3, 



h 3 



/ In 1 + e ^ s (^ s_c v /m i c2+ l?'l 2+I/sC P' a+ ' 7s ) 
</r 3 L 



dp 



for the semi-classical model, and by 



gf(u a ) = a; 



/l 3 



3 /3 s (/j s — c-y/ m|c 2 + |p| 2 +^ s cp-a+!7s) 



dp 



R 3 



(38) 



(39) 



for the classical model. Notice that in the classical model we have M s = /J" 1 ]^, while in the 
semi-classical model we have M s > f3~ 1 N s by the simple convexity inequality lnrr < —1+x, 
with "=" only for x — 1. Notice furthermore that, in either case, the map w i— > g s {w) is 
monotonic increasing. 

Lemma 4.1: Let the pair (u + ,U-) solve equations ffiSTty and (35), s = =F, under the 
asymptotic conditions with 7 given in satisfying $ZdJ. Then 

\Zs,tls,tN s N t - Z S M S > 0, (40) 
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and equality in (JQ) holds if and only if both u + and w_ are radially symmetric and de- 
creasing about the same point. 



Proof: We follow the general reasoning of [It], [XT], [12| . 

Since, by hypothesis, the pair (u + , uJ) solves the equations fl32| ) and (|33|), s = =F, under 
the asymptotic conditions (|34"|), then {U + ,UJ) satisfies (|29| ) and fl33|) , s = =F, under the 
asymptotic conditions fl3"T|). Therefore, as |x| — > oo, 

G s (C/ s )(x) = ^ / e /3a(M S -cV^ 2 +bl 2 +w») dp | x | -2/9. Et 7..* W0»)) ; (41) 

with #(x) = o(l). Also by hypothesis, ([3(]) is satisfied, so that from ([y]) we conclude 
that psJ^fisjNt £ 1. Then, by (|31]) again, and since U s G C 2 ' Q (hence, Ef, E C°° by 
bootstrapping), the level sets A| = > £ s } are compact, hence |A|| < 00. 

Let x 1 — > U*(\x\) denote the equi-measurable, radially symmetric, non- increasing rear- 
rangement of x t— > U s (x), centered at the origin, and denote by A|* = {x\ U* > £ s } the 
ball of radius r|, centered at the origin. By Sard's theorem the C°° regularity of the U s 

implies that the outward normal A to <9A| exists except at most for £- values in a set of 

measure zero, so that the ensuing manipulations involving A to <9A| are well defined £-a.e. 

First, recalling that G s > 0, we note that on <9A| we have — (A, VC/ S ) = |VC/ S |, by the 
Hopf lemma. Integration of this identity over 9A|, a trivial rewriting, and an application 
of the Cauchy-Schwarz inequality now gives the estimate 

(X,VU s )da= [ \VU s \ 2 -±—da>([ da )( f j^-rdr) , (42) 

9A| JdA° \vU s \ yJdA^ J \JdA^ \vU s \ J 

with equality holding if and only if \VU S \ is constant on 9A|. Noting that 

/ do- = |A|I, (43) 

JdA* 

and applying the classical isoperimetric inequality [0, we have 

|A«I > |0Af |, (44) 

with equality holding if and only if, up to translation, <9A| = c?A|*. By the co-area 
formula 

I f 1 

da = / — — da. (45) 
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Pulling these estimates together we have 



(X,VU s )da>\dAir[ I T^da) , (46) 



5A| \JdA s s 



with equality holding if and only if, (i), |VC/ S | is constant on c?A|, and (ii), <9A| = <9A|*, 
up to translation. This last remark implies in particular that we can restate fl46"[ ) as 0, 

- / (A, W s ) da>-[ d r U; da. (47) 

JdA° </<9A|* 



Next, using Green's theorem and (|29|), then rearrangement identity for s = t, then 
rearrangement inequality for s ^ t (in which case t = —s), noting that 7 S ,_ S > and 
recalling that w h- > G s (w) is increasing, we have 



(A, W s ) da =- AU S dx (48) 

,/A| 

= ^Y,tls,t I G t {U t )dx (49) 

</A| 

= Ak( 1s J G S (U;) dx + i at - s [G- s (U-s)dx) (50) 
\ Ja°* ja| / 

<^Hfts,tl G t {U* t )dx (51) 

J A** 

where equality in fl5"T|) can hold only if U t and U s share their level lines (up to the labelling) 
in A|, for our 7 is irreducible. 

Combining inequalities ([IT]) and fl5"T|), we arrive at the inequality 

- / 8 r U: da < ^Y.tls,t I G t {U*) dx, (52) 
7aA|* Ja°* 

where equality can hold if and only if each A| is a disk, with \ VU S \ constant on 9A|, and 
all the U s share their level lines (up to the labelling). Thus, in case of equality in (|52"D , 
from the first two conditions for equality it follows that the family of disks A^" and the 
family of disks A^ are separately concentric, while from the third condition for equality it 
then follows that the families of disks must be jointly concentric. 
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On the other hand, if at least one of the U s is not radially symmetric decreasing about 
any point, let S s be the image under U s of the (generally non-radial) set C M 2 which 
supports the non-radial parts of U s . Then S s has finite measure. Since equality in ( |5^ ) 
cannot hold for £ G S s , for £ G? we now conclude that we have strict inequality in fl52|), 

- 2wr\U?{r§ < 47r£ tTs , t / G t (U*) dx (53) 



for both s = =F. 
We now set 



and 



AA s (r)= / G s (U* s )dx, (54) 



M a (r) = / g s (U;)dx. (55) 

JB r (0) 

We have lim^oo J\f s (r) = N s and lim^oo M. s { r ) — M s , for 

f f(U:)dx= [ f(U s )dx, (56) 

where / stands for either g s or G s . By (|53"D, 

27rrC/;'(r) > -47r£ tTS)i M(r), (57) 

from which we conclude that 

rM" s {r) > M' s (r) - 2A/»£ t 7 M M(r), (58) 

with ">" valid for all r > for which £/*(r) G 2 s , while "=" holds for U*(r) G" E s . We 
now sum (|58"D w.r.t. s = =F, obtaining 

E/KW > E^.(r) - E M 7,,t(^(r)M(r))', (59) 

where we made use of the fact that 7 is real symmetric. Next we integrate (p3) from 
r = to r = 00, using integration by parts on the left-hand side. Since g s (U*) G L 1 (M 2 ) is 
radially decreasing, we have 

lim rM' s {r) = 0, (60) 
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thus we get the result 

lZ s ,tls,tN s N t - Z S M S > 0. (61) 



Now, if "=" holds in (|6"I), then all the level curves 9A| are circles with \WU S \ constant 
on c?A|; hence the circular level curves of each U s are concentric, and then U s (x) = 
U*(\x — Xq\) for some Xq. Moreover, since in case of "=" in ( KTf ) ( |5TD tells us that the 
two U s must share their level curves (with generally different level values, of course), we 
conclude that Xq = Xq , i.e. U + and LL are then radially symmetric and decreasing about 
the same center of symmetry, xq. 

On the other hand, if at least one of the U s is not radially symmetric and decreasing 
about any point, then the integration picks up all the strict inequalities from ( 6?, and 
">" holds in (gip. 

Finally, since rank (7) = 2, it follows that at least one u s is not radially symmetric 
and decreasing about any point if at least one U s is not. In the same vein, u + and u_ are 
radially symmetric and decreasing about the same center of symmetry, xq, whenever both 
U + and LL are. 

This proves Lemma 4.1. Q.E.D. 

4.3 Asymptotic control near infinity 

Standard harmonic analysis gives us: 

Proposition 4.2: Under the hypothesis stated in Lemma 4-1, each solution pair (u + ,U-) 
°f dHA satisfies the integral representation 

u s (x)-u s (0)= [ (ln r -^-ln T ^)G s (j2 tlS)t u t )(y)dy. (62) 

Jw.2 V \x-y\ \y\ J v } 

Corollary 4.3: By ffiiSJ we have 

V« s (x) = -2 / Gs (j2 tJs>tUt ) (y)dy. (63) 

J R 2 \x — y\ \ / 



With the help of (B3) and (^3) we obtain asymptotic control over the r.h.s. of (|3"2"D, 
expressed in terms of the U s . 

Lemma 4.4: Under the hypotheses of Lemma 4-U there exists an ro(U s ) > 0, a constant 
C s > 0, and a monotonic decreasing h s (\x\) > satisfying 

lim \x\- hsilxl) = 0, (64) 

|a;|^oo 
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such that for s = =F we have, for \x\ > r^, 

G s (U s )^)<Cs\x\- 2 - h ^ xl) . (65) 

Furthermore, for at least one s, we have h s (\x\) > e s > for \x\ > r . 

Proof: The bound (|65|), with /t s (|x|) = 0(1) monotonic decreasing, follows directly from 
(flip and fl30|) . Furthermore, by (fill) and ( |30"D we can find h such that |x|~ 1_h ^' :r ^d|x| < 



oo, but this is impossible if |x| _/ls( -' a ''' ) -f* 0; hence, (B3) follows. 

This still allows /i s (|x|) = o(l) for both s, but by Lemma 4.1 and the fact that M s > 
fij x N a (see the definition of the g s above), we find, after multiplying (|6lf) by —2 and 
re-grouping terms, that 

J2 s P7 lN s (2 - AEtT./i) < 0- ( 66 ) 
Thus, if for one of the s we have /i s (|x|) = o(l), say for s = +, then by (fill) we have 
/3 + 52*7+,*-^* = 1, so that ( j66|) gives us right away 



/3-E t 7-, t AT, > 2 + (67) 



By symmetry, the analog conclusion holds for fl + ^2 t j +tt Nt if /i_(|x|) = o(l). Hence, 
/i s (|x|) = o(l) for at most one of the s. 
This proves Lemma 4.4. Q.E.D. 

Corollary 4.5: For at least one s, we have \J R2 ln\y\ G s (U s )(y)dy\ < oo, so that for this 
s, we have 

lim (u s (x) + 2N s \n\x\) = u s (0) + 2 [ \n\y\G s (U s )(y)dy. (68) 
|ai|->oo v / y R 2 



We proceed with gradient estimates. 
Lemma 4.6: Under the hypothesis stated in Lemma 4-1, each solution pair (iz + , uJ) of 



(pq), (\33j), ( Ww satisfies the gradient estimates 



lim sup \x\\Vu s \ < 2N S . (69) 

[a;[— »oo 



Proof: By Corollary 4.3, we have 



|V//,(,-)| <2 / G f s) dy. (70) 
\x — y\ 
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After multiplying flTO|) by \x\, a simple rewriting of the r.h.s. gives 

\x\\Vu s (x)\<2 [ G s {U s ){y)dy + 2 [ ( - l) G s {U s )(y)dy. (71) 
Jr2 J R 2 \\x-y\ J 



By ( p0|) the first integral on the r.h.s. of (|7T| ) equals N s . By the triangle inequality, the 
second integral on the r.h.s. of ([71]) is bounded in absolute value by 

2 f r \y\ Gs (U s )(y)dy. (72) 

Jr 2 f — y\ 

We now show that 

lim / G s (U s )(y)dy = 0, (73) 

from which the lemma follows. 

We split the domain of integration in (f72"|) as follows: 1R 2 = fiiUf^Uf^, with Qi = {y \ 
\y\ < \x\/2}, n 2 = {y | \x\/2 < \y\ < 2\x\}, and fi 3 = {v I 1 2/ 1 > 2|z|}. If G a (U a )(y) < 
C\y\~ 2 ~ e , with < e < 1, then the estimates are precisely the same as in |TI|, section 2, 
with exp replaced by G s ; this is the case for at least one of the s. It remains to provide 
estimates when ft.«(|a;|) = o(l) for one of the s. 

To estimate the contribution from Oi when G s (U s )(y) < C\y\~ 2 ~ ha ^ v » with /i s (|j/|) = 
o(l), we note that 

/ ^-G 8 {U s ){y)dy<^- f \y\G s (U s )(y)dy < ^ C MC) dC (74) 
1^ 2/1 Fl Fl </o 

As for the right hand side of (0), L'Hopital's Rule gives us 

(~« f\ x \ 

lim — / C _/ls(0 dC = C lim |x|"MN) = 0, (75) 

\x\^oo \x\ Jq |a;|— +oo 

the last step by Lemma 4.4. Hence, the Lh.s.(fT4T) vanishes as \x\ — > oo. 

Similarly, the contribution from f2 2 is estimated by using again that G s (U s )(y) < 
C\y\~ 2 ~ hs ^ y \\ so that for \x\ large enough we have the bound 

/ r^GsiU^dy < C f pi < cr^i-MW). ( 76 ) 
Clearly r.h.s. ([781)— > as |x| — > oo, by the same reasoning as for f^. 
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Finally, the contribution from is dominated by 



I V I 



;G s {U s ){y)dy<C G s (U s )(y)dy, 
o 3 \ x ~~ y\ J\y\>z\x\ 

which vanishes as |x| — > 00, by hypothesis (|30"1). 
This concludes the proof of Lemma 4.6. Q.E.D. 

Lemma 4.7: Under the hypotheses of Lemma 4-1, we have, uniformly in x, 

lim (x, Vu s ) = -2N S . 

\x\— >oo 



(77) 



(78) 



Proof: Let x = x/\x\ and y = y/\y\, with \x\ = \y\. Now fix x G § . By (|34|), we have 

u s (tx) 



lim 

t^oo In r 



-2iV„ 



(79) 



Thus, by L'Hopital's Rule, 



lim t—u s (tx) = lim (1, V« s ) = —2N S 

r^oo dr \x\— >00 



(80) 



for x = \x\x. It remains to establish uniformity of fl80[ ). To this extent, we show that there 
exist R and 5, such that, if \x\ > R and \x — y\ < 5, then 



\(x,Vu s {x)) - (y,Vu s {y)) \ < e. 
We first show that for \x\ > R and \x — y\ < |x|/10, we have, 

\x\\Vu s (x) - Vu s {y)\ < C\x -y\ + C\x\- hs ^\ 

By Corollary 4.3, 

1 x — z y — z 



\Vu s (x)-Vu s (y)\<2 / G s (U s )(z 



\x — z\ 



\y - z \ 



dz. 



(81) 



(82) 



(83) 



We break up the domain of integration in the above integral exactly as in the proof of 
Lemma 4.6. (Notice the integration variable is now z.) The integration over Qi is estimated 



exactly as in section 2 of [10] to get 



G s (U s )(z) Kj±dz < cKJ± 



X 



\X\ 



(84) 
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The integral over is dominated by 

/ G s {U s ){z) l-r-^—r + r^—^) dz ^ C\x\- l - h ^\ (85) 

J\z\~\x\ \\ x ~ A \y~ z \J 

The final estimate above was identical to that made in the proof of Lemma 4.6. Use was 
made of G s (U s )(z) < C\x\~ 2 ~ h ^ on 0,2, which holds by Lemma 4.4. The contribution 
from 3 is estimated once again exactly as in section 2 of pi to be dominated by 



C\x - y\ ! Ga ™ z) dz < C'^^ [ G s (U s )dx < C"^^, (86) 

J\z\>2\x\ \ z \ \ x \ JK? Fl 

where the last step follows by ([ft]). By these estimates, 

\(x,Vu s (x)) - (y,Vu s (y))\< \x\\x - y\\Vu s (x)\ + \y\\Vu s (x) - Vu s (y)\ (87) 

< \x\\Vu s (x)\\x -y\ + \x -y\ + C\x\- h 'W\ (88) 

By Lemma 4.6, the last expression above is at most C'5 + C\x\~ hs ^ x ^ . Thus our claim 
( [81]) follows now from Lemma 4.4 for suitably large R and small 5. Since S 1 is compact, 
uniformity of the limit in Lemma 4.7 now follows. Q.E.D. 

Lemmata 4.6 and 4.7 imply 
Corollary 4.8: Under the hypotheses expressed in Lemma 4-1, we have, uniformly in x, 

lim \x\\Vu s \ = 2N S . (89) 



Proof: Follows essentially verbatim fllQfl , proof of Corollary 2.2, with exp replaced by G s . 
Q.E.D. 

Let = {x\ u s (x) > £}, where £ <C —1. By (|34|) it follows that if x G cft\, then 
|x| > R(c) with R(c) large. For such x, it follows from Corollary 4.8 that Vu s ^ 0. Since 
u G Cf£ we easily see that therefore dQ% G C 2,a . Thus the unit outward normal &{x) to 
dfl^ exists at all x G dQ^ for £ sufficiently negative. 

Lemma 4.9. Let uj(x) be the unit outward normal to dfl^ at x, and let x = x/\x\. We 
have, uniformly in x, 

lim (x,cu) = 1. (90) 



Proof: Identical to |10| , proof of Lemma 2.8. Q.E.D. 



Remark: Lemma 4.9 implies that asymptotically for large x the dQ^ become concentric 
circles. 

We are now in a position to prove our Theorem 3.1. 
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4.4 Proof of the Virial Theorem 

Proposition 4.10: (Rellich-Pokhozaev identity.) Under the hypotheses expressed in The- 
orem 3.1, any solution pair (u+,U-), of flSlfy, ftSDj), ffift ) satisfies the Rellich-Pokhozaev 
identity 

\j2 s ^s,tN s N t - Z s Ms = 0. (91) 

Remark: The Rellich-Pokhozaev identity fl9~T| ) is identical to the identity expressed in the 
Virial Theorem 3.1. 



Proof of Proposition 4-10: For (u + ,m_) a solution pair of (|32]), (p3|), (|34|), we have the 
partial differential identity 

div«ac, Vut) Vu s ) = (Vu„ (1 + (x, V»V«t> - 4tt(x, Vu t )G s (E t 7 s ,^)- (92) 

We will multiply fl92|) by 7 si , sum over s and t, integrate over _B^, use some partial 
integrations, then take the limit R — > oo. 

We evaluate first the left-hand side of (|92|). Green's theorem gives us 

/ div((x, Vu t )Vu s )dx = / Vu t ) (x, Vw s ) da. (93) 

J B R JdB R 

Taking the limit R — > oo, using (|78f) , we get 



hence 



lim / Vu t )(x, Vu s ) da = 8irN s N t} (94) 

JdB R 



lim E st 7s,t / div((x,VM t )Vu s )dx = 87r^ i7sit Ar s Ar t . (95) 
On the other hand, the last term in the right-hand side of (^2j) gives us 

J2 s ,ffs,t ( x > VM t )G s (X] t 7s,tMt)dx = Y,s / ( x > ^9s(J2ffs,tU t ))dx 

J B H J B R 

/ M 5 , s (E i 7 s ,tMt)do- - 2]T s / ^ s (Xl t 7« ) t«t)da ? - (96) 
JdB R Jb r 
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We now take the limit R — > oo in (|96|). As for the surface integrals, we note that by 
Lemma 4.4 we have g s {U s )(x) ~ CG s (U s )(x) as \x\ — > oo, so that once again by Lemma 



4.4, we have 



As for the volume integrals, we get 



lim / \x\ g s ^2fis,tUt) da = 0. (97) 

l -™ JdB R 



lim / g s (J2 t ls,tU t ) dx = M s . (98) 
R ->°° Jb r 



Turning now to the first term in the right-hand side of (|92|), we use the symmetry of 7, an 
integration by parts and (|32|) , to get 

E s ffs,t I (Vu s ,(l+(x,V))Vu t )dx = V st 7M / \x\(Vu s ,Vu t ) da (99) 

JB R 1 JdB R 

= \llst1s,t [ (\x\- 1 (x,Vu s )(x,Vu t ) + \x\(V T u s ,V T u t ))da (100) 

1 ' JdB R V ' 

where Vt denotes tangential derivative. By Lemma 4.7 and Corollary 4.8, we have 

| x | 2 |V T u s | 2 = |x| 2 |Vu s | 2 - (x, Vu s ) 2 -> 0, (101) 
uniformly as \x\ — > 00. Thus as R — > 00, 

|z|<V T « a , V T w t }do -> (102) 



and therefore 



lim E s t7»,t / (Vu s , (1 + (x, V))VtZt) dx = 4vr£ s t 7.,t#-#t- (103) 

Pulling all limit results together, we obtain Proposition 4.10. Q.E.D. 

Remark. The proof of the virial theorem extends to situations when 7 does not have full 
rank, hence to more-than-two species beams. 
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4.5 Concluding the proof of the Symmetry Theorem 

By Lemma 4.1, and by Proposition 4.10, the solutions u s of (|32|) , (|33"D , (|34"1), have to be 
radially symmetric and decreasing about a common center x$. Since the coupling matrix 
7 is invertible, the same conclusion holds for the solutions U s of (|2~9"D, ([301), fl3~lf ). The proof 
is complete. Q.E.D. 
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